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ABSTRACT 

In a sample of 54 galaxy clusters (0.04 < z < 0.15) containing 3551 early- type galaxies suitable 
for study, we identify those with tidal features both interactively and automatically. We find that 
~ 3% have tidal features that can be detected with data that reaches a 3<r sensitivity limit of 26.5 
mag arcsec -2 . Regardless of the method used to classify tidal features, or the fidelity imposed on 
such classifications, we find a deficit of tidally disturbed galaxies with decreasing cluster centric radius 
that is most pronounced inside of ~ 0.5i?200- We cannot distinguish whether the trend arises from 
an increasing likelihood of recent mergers with increasing clustercentric radius or a decrease in the 
lifetime of tidal features with decreasing clustercentric radius. We find no evidence for a relationship 
between local density and the incidence of tidal features, but our local density measure has large 
uncertainties. We find interesting behavior in the rate of tidal features among cluster early- types as 
a function of clustercentric radius and expect such results to provide constraints on the effect of the 
cluster environment on the structure of galaxy halos, the build-up of the red sequence of galaxies, and 
the origin of the intracluster stellar population. 



1. INTRODUCTION 

Galaxies and clusters of galaxies assemble hierarchi- 
cally, through a long progression of mergers and ac- 
cretion. As one considers galaxy-scale dark matter 
halos, baryonic physics becomes more integral and so 
the analytic formalism develop ed to calculate the rate 
of dark matter halo evolution ([Press fc Schechterlfl974t 
iBond et all 119911: lLac ev fe ColeTTT993^ which is broadly 
confirmed by numerica l simulations (jLacev fc C ole 1994; 
iSheth fc Tormenl Il999[ ). becomes less prescriptive. To 
remedy this shortcoming, galaxy formation is now stud- 
ied usi ng hydrodynamic nu merical simulations (e.g., 
iSpringellBOOa iHoffman et al.ll2010D . although the details 
of the assembly process are sensitive both to baryonic 
physics below scales accessible in the simulations ("sub- 
grid" physics) and to the complex history that occurs on 
cosmo logical scales. Tests of such models include com- 
parisons to the properties of galaxy clusters, their con- 
stituent galaxies , and the generation of intracluster stars 
(jPuchwein et al.ll2010D . 

The rate of accretion events and mergers is a funda- 
mental quantity that hierarchical growth models must 
be able to predict accurately. A proper compari- 
son to an empi rical rate measurement is elusive (see 
lLotz et all 120111 for a description of the current state 
of the field). Efforts to measure this rate, particu- 
larly as a function of redshift, have focused either on 
meas uring close pairs of galaxies (a few examples in- 
clude Carlberg, Pritchet, & Infante 1994; Patton et al. 
19971: iLe Fevre et al.l 120001: iKartaltepe et al.l 120071 : 



de Ravel et al .1120091 ), or identifying galaxies that appear 
morphologically to be the result of on going or recent 
merger events (a few examples include Abra ham et al.l 



1996; lLotz et al.l 2008: J ogee et al.l 120091 : iBridge et all 
2010; iMiskolczi et al.l |2011[) and results often disagree 
among studies (jLotz et al.ll2011h . These measurements 
are difficult, even with the superb angular resolution 
of the Hubble Space Telescope, for high redshift galax- 
ies. However, renewed focus on identifying tidal fea- 
tures at low redshift has yielded some striking suc- 
cesses botli_jntheJncreased incidence of detected fea- 
tures (|Tal et al. 2009) a nd the ir scale and morphologies 
(jMartmez-Dekado et al.ll2010D . 

We aim to make a measurement of the major merger 
rate at low redshift, but as a function of environment. A 
complication in practice is that the definition of a major, 
vs a minor, merger is typically made on the basis of the 
mass ratio between the two galaxies. Such a definition 
is difficult to apply empirically, particularly in environ- 
ments that may have affected the individual galaxies and 
in cases where the two initial galaxies have merged into 
a single remnant. We define a major merger here as 
one that results in the large-scale tidal features we are 
searching for around luminous, early-type (R 1 / 4 profile) 
galaxies. Although this approach evidently complicates 
comparisons to models, it should not impact our primary 
goal of searching for patterns in the occurrence of tidal 
fe atures as a function of environment. 

iDubinski, Mihos, fc Hernquistl (|1999f ) showed how the 
dark matter halo structure of the merging galaxies can 
significantly alter the appearance of a major merger, 
thereby demonstrating that measurements of major 
merger signatures, even at a single redshift, could be a 
sensitive model diagnostic. Beyond the purpose of under- 
standing hierarchical accretion better, we need to study 
the merger rate in clusters and in the surrounding en- 
vironments because these events contribute stars to the 
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ubiquitous intracluster stellar population (IZibetti et al.l 
2005; Gonzalez, Zabludofl, & Zaritsky 2005}, which is an 
important component i n estimates of the baryon budget 
of gr o ups and clusters (E onzalez, Zaritskv, & Zabl udofl] 
2007; Giodini et al. 2009; McGee & Balogh 2010), and 



intracluster medium ( 


Zaritskv. Gonzalez. & Zabfudofll 


2004; Sivanandam et al. 


2009). Furthermore, extending 



any study of galaxies within clusters out to the surround- 
ing environs is critical because differences in their star 
formation histories, which could be related to interac- 
tions, be£in_to_jppear at cluster centric radii of several 
Mpc (jLewis et al.l 120021 : iGomez et al.ll2003[ ). 

The principal observational challenge in identifying 
tidal features is that such features are of low surface 
brightness and have irregular morphologies. A vari- 
ety of techniques have been devel oped to identify merg- 
ing galaxies and mergers (e.g., Conselice etldl 120031 : 
lLotz et aTl [2008), but the problem takes on different 
forms depending on the redshift of the galaxies and the 
depth and uniformity, or "fla tness", of the data. For 
galaxies in the local universe, IColbert et all (2000) em- 
ployed unsharp masking (subtracting a smoothed image 
of the galaxy from itself), galaxy model division (dividing 
the galaxy image by a best-fit model galaxy made from 
nested elliptical isophotes), and color mapping (division 
of the galaxy image in one color band by the galaxy im- 
age in another color band) to aid in the visual detection 
of tidal features and found that 41% (9 of 22) of analyzed 
isolated galaxies show tidal features, while only 8% (1 of 
12) of group galaxies had such features. While this re- 
sult is certainly suggestive of lower rates of tidal features 
in group environments than in the field, Poisson statis- 
tics limits the sign i ficance of the discrepancy to less than 
2cr. IVan Dokkuml ([20051 ) divided his galaxy images by 
a model fit and measured the mean absolute deviation 
of the residuals to define a quantitative threshold for the 
identification of tidally disturbed galaxies. He found that 
53% of his sample of 126 red, field galaxies show tidal fea- 
tures, while 71% of the bulge-dominated early- type sub- 
set of 86 galaxies appear to be tidally disturbed. Using 
a similar met hod for a sample of nearby, luminous, ellip- 
tical galaxies, iTal et al.l J2009) found that 50% (5 of 10) 
of cluster galaxies, 62% (13 of 18) of poor-group galax- 
ies, and 76% (16 of 21) of isolat ed galaxies ha v e tida l 
features. Similar to the results of iColbert et al.l ((2000), 
these also suggest that tidal features are less common 
in cluster environments than in the field, but again the 
significance of the resu lt is limited by the sample size. 
iJanowiecki et al.l (|201Q[ ) quantify the extent of tidal sub- 
structure in five elliptical galaxies in the Virgo cluster by 
measuring the luminosity of the substructure and find 
no obvious correlation between clustercentric radius and 
the amount of substructure. Larger samples are needed 
to resolve this issue. 

Large r samples studied thus far are from field sur- 
veys. Bridg eet al.l (j2010h visually identify tidal tails 
in a sample of 27,000 galaxies over 2 square degrees of 
the Canada-France-Hawaii Telescope Legacy Deep Sur- 
vey (CFHTLS-Deep) to find that the merger fraction of 
galaxies evolves from 4% at z ~ 0.3 to 19% at z ~ 1. 
Miskolcz i et al.l (|2011f ) study a sample of 474 galaxies se- 
lected from the SDSS DR7 archive to find that at least 
6% of the galaxies have distinct tidal streams and a to- 



tal of 19% show faint features. No large samples have 
included substantial numbers of galaxies in dense envi- 
ronments. 

Deep, uniform imaging is a prerequisite for any such 
study. The incidence of identified mergers among field 
galaxies incr eases from < 10% loc ally in relatively shal- 
low imaging (Mis kolczi et al.ll2Qllh to many tens of per- 
cent in deeper imaging (Martmez-Delga do et al.l 2010; 
ITal et al.ll2009l ). For this and other reasons, it is nei- 
ther surprising nor incorrect that the apparent total rates 
from different studies vary widely. The key is therefore to 
compare within a given study rather than across studies. 

We identify tidally disturbed galaxies in the largest ex- 
isting sample of deep, wide-field images of nearby galaxy 
clusters (0.04 < z < 0.15). These come from the Multi- 
Epoc h Nearby Cluster Survey (MENeaCS; iSand et al.l 
2012), and this study is part of a series using those data 
to address a series of science questions including so far th e 
incidence of intracluster super novae (ISand et al. 2"0lTh . 
the rate of SNe la in clusters (Sand et al. 2012), t he evo - 
lution of the dwarf-to-giant ratio dBildfell et al.lj20TT 
and the rate of SNe IPs in clusters (jGraham et al.ll201^ 
The larger sample size enables us to examine, for the first 
time, relative rates of tidal features as a function of dif- 
ferent measures of the global environmental (projected 
clustercentric radius) and local environment (projected 
local galaxy density). In §2 we briefly present the data 
we use and discuss our analysis of these data, including 
the selection of candidate cluster galaxies and the calcu- 
lation of the tidal parameter. In §3 we present our results 
and discuss whether any significant trends are identified. 
We present our conclusions in §4. We adopt Ho = 71 km 
s _1 Mpc -1 , Qm — 0.27, and Q\ = 0.73 for conversions 
to physical units. 

2. DATA AND ANALYSIS 

2.1. The Cluster Sample 

The galaxy cluster images used for this study are 
the deep stack images obtained by the Multi-Epoch 
Nearby Cluster Survey of X-ray selected clusters (ME- 
NeaCS; [SanJ^F2]|20T2|) with MegaCam on th e Canada- 
France-Hawaii Telescope (jBoulade et al.ll2003h . The im- 
age stacks typically contain 20 to 30 120-second expo- 
sures in the g' and r' filters for each cluster, observed at a 
monthly cadence t o pur sue the supernovae science goals 
(jSand et al.l l201ll l2012f ). The analysis is done on the 
deeper r' data. The depth of the r' image stacks for point 
source detection corresponds to limiting absolute magni- 
tudes of —11 to —13 at the distance of the clusters. The 
3cr limiting apparent surface brightness, derived from the 
sky RMS, of 1 arcsec 2 patches in the deep stack images 
is approximately 26.5 mag arcsec -2 . The faintest tidal 
features that we visually identify are approximately this 
magnitude (see Fig. [71 and [8]). MegaCam features a large 
field of view (~1 degree 2 ) that enables us to search for 
tidal features from the central regions out to 1 — 2i?200 
and study the radial and environmental dependence of 
the merger rate. Details on the clust er parameters and 
imag e reduction process are given in San d et aTl (2011, 
120121 ) . The list of the 54 clusters analyzed is given in Ta- 
ble [j The cluster redshifts range from 0.04 < z < 0.15, 
with the redshifts taken from the NASA Extragalactic 
Database. 
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TABLE 1 

Galaxy Clusters 



Cluster 


Redshift 


^X,0.1 

(10 44 


-2AkeV 

ergs/s) 


M 20 o a 
(M ) 




#200 b 
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Visual Sample 


Ngal 
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5.12 


7.28 
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72 


2 


1 
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99 


11 


2 
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7.44 
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8 


2 
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32 





3 
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93 


15 


4 
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38 


2 


1 
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3.85 
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71 


3 


2 
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Abelll3ol 


0.117 
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57 


2 


1 
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0.143 
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72 


5 
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44 


5 


2 
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3.79 
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30 


1 


2 
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62 
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26 


2 


2 
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17.44 


16.57 


2380 


75 


4 


2 
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2.55 


4.38 


1530 


76 


2 





Abell2050 
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4.54 
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86 


6 


3 
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1670 


54 


1 


1 
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55 


1 
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1 


Abell21 


0.095 




2.64 


4.51 


1540 


43 


2 


4 


Abell2142 


0.091 




21.24 


19.56 


2510 


130 


1 


2 


Abell2dl9 


0.056 




15.78 


15.61 


2350 


24 


3 


2 


a i ii r> j r>r\ 

Abell24U9 


0.148 




7.57 


9.69 


1950 


88 








Abell242U 


0.085 




4.64 


6.67 


1760 


64 





3 


Abell242o 


0.098 




4.96 


7.04 


1780 


88 





1 


Al no /i /in 

Abell2440 


0.091 




3.36 


5.33 


1630 


63 








A 1_ „ 1 1 O /I /I O 

Abell2443 


0.108 




3.22 


5.22 


1600 


46 





1 


Al no/inr 

Abell2495 


0.078 




2.74 
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Note. — Basic properties of clusters analyzed along with the total number of galaxies and number in the 
visual and cleaned-auto samples of disturbed galaxies for each cluster. 
a M200 derived from L x ,o.i-2AkeV using relation found by Reiprich & Bohringer (2002) 
b #200 found from correspondence with M200 



2.2. Selecting Elliptical Cluster Galaxies 

We use SExtr actor, a source detection algorithm 
(Bert in fc Arnoutil 19961 ), in two-image mode to measure 
the magnitudes and colors of all galaxies in the images. 
The galaxy clusters are relatively nearby (z < 0.15) and 
we can only detect tidal features in galaxies of relatively 
large angular extent (many arcsec), so we select for fur- 
ther study only sources that have a minimum area of 60 
pixels that are each at least 6a above the background 
and have a stellarity parameter of less than 0.05. Using 
that output, we create color magnitude diagrams for each 
cluster using the absolute magnitudes derived from SEx- 



tractor's AUTO-MAG parameter , which uses an aper- 
ture radius of 2.5rx ron (|Kronlll9"80l ). and the color using 
the APER_MAG parameter, which finds the magnitude 
within a small aperture we set to 1.9" in diameter in or- 
der to maintain high S/N and limit contributions from 
neighboring sources. 

To identify a sample of galaxies with a high proba- 
bility of being at the cluster redshift and to avoid star- 
forming, late-type galaxies that can have significant in- 
ternal structure, we select red sequence galaxies for fur- 
ther study. To do this, we first visually examine the 
location of galaxies in the color-magnitude diagram that 
are within 1.5 arcmin of the cluster center. Using these 
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■ ■ All Galaxies 

♦ ♦ Central Galaxies 

I I Selection Window 



Fig. 1. — Sample color magnitude diagram (Abell 21). The black 
dots represent all galaxies meeting our SExtractor selection param- 
eters. The red diamonds are galaxies located within 1.5' of the 
cluster center used to originally define the location of the cluster 
red sequence. Galaxies within the shaded region are brighter than 
our minimum magnitude and are identified as red sequence mem- 
bers of the cluster. These are the galaxies in which we search for 
tidal features. 



galaxies, we define the color normalization of the red 
sequence and adopt an empirically-determined slope of 
—0.033. The red sequen ce slopes of clusters at a given z 
show little if any scatter ([Ellis et al.|[l9 97: Stanf ord et al.l 
19981 : IGilbank et al" 2008). We select galaxies with colors 
within ±0.13 of the mean relation for further study. This 
estimate of the color s catter is probably so mewhat con- 
servative (for example, Bildfel fet al.l (j2Q12h select galax- 
ies within ±0.2 of the red sequence), but it minimizes 
contamination, which for this study is more important 
than completeness. A tight color selection could result 
in a bias if interacting galaxies are bluer than the mean, 
but (as we find later) they are not. We do not explicitly 
apply k-corrections or evolution corrections, but they are 
inherent in our red-sequence color selection. We repeat 
the procedure of visually identifying the optimal zero 
point for each cluster, keeping the slope and thickness 
of the color-magnitude band constant. We then apply 
these cuts to the entire image, keeping only galaxies with 
M r < — 20 (where M r is simply the sum of the appar- 
ent magnitude and the distance modulus) at the cluster 
redshift (see Fig. [1] for a sample color-magnitude dia- 
gram and galaxy selection) to ensure that the candidate 
galaxies are sufficiently luminous for our tidal analysis. 
We will discuss estimates of the interloper fraction in 
§3.2. 

We focus our study on elliptical galaxies because their 
surface brightness distributions can be modeled more 
easily than those of spirals and because disk galaxies are 
believed to not have been involved in major mergers since 
moderate redshifts unless both progenitors were gas rich, 
a conditio n that will be rare in t hese high density envi- 
ronments ([Robertson et al.ll2006[ ). To limit our sample 
to elliptical galaxies, we reject galaxies with ellipticity 
greater than 0.5 and we discriminate based on the shape 
of the surface brightness profile. To do the profile dis- 
crimination, we select the radial region of the galaxy's 
surface brightness profile in which the intensity is less 
than 10% of the galaxy's peak intensity (to avoid the 
exponential bulge of spirals) and greater than la above 
the background sky level. We reject galaxies that have 




20 40 60 80 1 2 3 

SMA (in pixels) SMA 1 / 4 (in pixels) 

Fig. 2. — Sample luminosity profile fit used to identify a spi- 
ral galaxy. The linear least-squares fit of the In (intensity) versus 
semi-major axis (left panel) has a smaller % 2 profile than the least- 
squares fit of the ln(intensity) versus semi-major axis 1 / 4 . The por- 
tions of the luminosity profile used for the best fits are shown in 
red, the portions excluded are shown in black. The blue lines are 
the best fit linear regression lines. 



a lower \ 2 f° r the weighted least squares linear fit of 
ln(intensity) vs. radius than of ln(intensity) vs r 1 / 4 (Fig. 
[2]). Out of the initial sample of 11,904 red sequence galax- 
ies, 4526 were flagged as disks due to their surface bright- 
ness profiles and 2322 galaxies were flagged as disks due 
to their ellipticity, with a union of the two criteria result- 
ing in a total of 5486 galaxies being dropped from the 
initial sample. We found with visual examination that 
only ~ 3% of galaxies classified as ellipticals had obvi- 
ous spiral structure (these are later discarded from the 
final sample). It is more difficult to estimate the fraction 
of galaxies identified as disks that are actually ellipticals 
due to the similar morphology of E and SO galaxies, but 
the vast majority of the galaxies identified as disks did 
visually appear to be correctly flagged. There were some 
instances (< 1% of all galaxies) where obvious mergers 
were classified as disks, but these tended to be galaxies 
that still had an obvious disk component and/or lacked 
a clear nucleus. 

For the subsequent analysis, we extract image regions 
of width corresponding to seven times each galaxy's di- 
ameter (as determined by SExtractor). The sizes of these 
images range from 46 arcmin 2 for the largest central 
galaxies to 500 arcsec 2 for the smallest galaxies that meet 
our thresholds, with the median size being 1 arcmin 2 . 
The angular sizes translate into 850x850 kpc to 50x50 
kpc with the median area being 105x105 kpc. The sizes 
of the images were chosen to ensure that any possible ex- 
tended tidal tails would be included in the analysis and 
that in almost all cases there is sufficient area to accu- 
rately define a background sky level. Visual examination 
did not reveal any cases where tidal features began be- 
yond the analyzed image region. Though there are a 
few cases where large tidal features continue beyond the 
analyzed image region, these galaxies are still correctly 
identified as tidally-disturbed by both the visual and au- 
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Fig. 3. — Flow chart showing process by which samples of tidally- 
disturbed galaxies are identified. 

tomated selections. 

2.3. Calculating the Tidal Parameter 

In studies of mergers and merger remnants there is 
always a tension between visual and automated iden- 
tification. Visual approaches have the advantage that 
they harness the tremendous pattern recognition skills 
that the human eye has evolved, but suffer from non- 
uniformity among different practitioners, and poorly de- 
fined criteria. Automated approaches are repeatable and 
well-defined, the latter being particularly important for 
comparisons to models, but are likely to be less powerful 
in selecting varying morphologies within this ill-defined 
class of objects. We resolve this conflict by applying both 
approaches. 

To objectively determine whether a galaxy is tidally- 
disturbed we fo llow the general procedure presented by 
iTal et all ()2QQ9h (shown in Fig. [3]) to calculate a tidal 
parameter from the average residual of the galaxy image 
divided by its best fit model. First a background gra- 
dient is fit and subtracted from each galaxy image us- 
ing the BACKGROUND generated by SExtractor with 
BACK_SIZE = 64 and BACK_FILTERSIZE = 3, in or- 
der to minimize large scale gradients across the galaxy 
images caused by intracluster light or nearby bright 
galaxies. The BACKGROUND map generated by SEx- 
tractor is a bicubic spline interpolation over the sigma- 
clipped pixels in the area of the BACK_SIZE parameter 
median-filtered over the BACK_FILTERSIZE number 
of BACK_SIZE areas. The BACK_SIZE parameter is 
chosen for our pixel scale (0.187" /pixel) to be larger than 
the scales of typical tidal features. 

The tidal parameter analysis is very sensitive to 
unmasked sources, which are easily confused with 
tidal matter. Pixels with known objects are masked 
using the SExtractor object catalog. To opti- 
mize masking of background sources while minimiz- 
ing masking of real tidal features, we create our 
masks from the union of two different masking cat- 
alogs. One catalog targets large galaxies using 
DETECTMIN_AREA = 5, DETECT_THRESH = 1.8, 
and DEBLEND_MINCONT = 0.001 to only mask very 
significant sources well removed from the subject galaxy. 
These parameters require DETECTMIN_AREA num- 
ber of neighboring pixels to be DETECT_THRESH 
sigma above the background in order to be detected as 



an object. The DEBLEND_MINCONT parameter de- 
termines if a cluster of pixels above the detection thresh- 
old should be divided into multiple sources, and has been 
tuned here to differentiate only very distinct clumps as 
separate objects. A second, more aggressive catalog tar- 
geting stars is generated by only identifying objects with 
CLASS_STAR less than 0.05 with DETECT_THRESH 
= 1.4, and DEBLEND_MINCONT = 0.0001 (a lower de- 
blending contrast to more readily identify stars blended 
with the target galaxy). We also create a back- 
ground sky frame to be used in a noise correction term 
(described below) using a third, very agressive cat- 
alog generated with DETECT_THRESH = 0.8, and 
DEBLEND_MINCONT = 0.001 that masks all sources 
(including the target galaxy and associated tidal fea- 
tures). Unmasked pixels are randomly chosen to replace 
the masked pixels in the creation of a background im- 
age in order to preserve the noise characteristics of the 
original image. 

We create a best-fit model for each galaxy using the 
IRAF tasks ellipse and bmodel. We fix the isophotal 
center but allow for variable position angle and elliptic- 
ity as a function of semi-major axis to generate best-fit 
isophotes out to 80% of the postage stamp image dimen- 
sions, and then create the corresponding model galaxy. 
The model image covers the entire postage stamp area, 
with pixel values set to zero beyond the area covered 
by the best fit model. Pixels whose model value ex- 
ceeds 0.01 times the peak value at the galaxy center 
are masked so that the analysis focuses on any extended 
tidal features rather than residuals near the galaxy nu- 
cleus that may arise from slight misplacement of the 
galaxy center in the model. The masked galaxy im- 
age is then divided by the model image. The resid- 
ual image is median filtered with a 5x5 pixel kernel 
in order to facilitate visual identification of tidal fea- 
tures. We further supplement our pixel mask by run- 
ning SExtractor on the smoothed residual image with 
DETECTMIN_ AREA = 5, DETECT_ THRESH = 5, 
and DEBLEND_MINCONT = 0.0001 primarily in or- 
der to mask faint galaxies close to the primary galaxy. 
The final mask for a typical galaxy includes 5-20% of 
the pixels in the area over which the model is computed. 
We discard all galaxies that have over 50% of the im- 
age area masked. This should not bias the results unless 
real tidal features are being masked (which we visually 
verified is very uncommon, with only a few cases out of 
several thousand galaxies). We also confirm with a con- 
trol sample that the masking is not responsible for the 
radial trend we discuss below (see §3.2). With the en- 
hanced mask we repeat the process of creating a model 
galaxy and smoothed residual image. 

We examine all galaxy and residual images for quality 
and calculate the tidal parameter. We visually inspect 
and discard the galaxy and corresponding residual im- 
ages that are contaminated with unmasked diffraction 
spikes from bright stars. We also flag galaxies that have 
visible tidal structure. This latter sample is referred to 
as our visual sample. We then generate an auto-detected 
sample using a tidal parameter that is calculated by find- 
ing the mean of the absolute value of the final residual 
image: 
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Fig. 4. — Example of galaxy in Abell 795 discarded interactively 
due to unmasked diffraction spike. Left column: galaxy image dis- 
played with 15" scale bar. Right column: residual of galaxy image 
divided by model with galaxy core and neighboring external ob- 
jects masked. The tidal parameter used to select tidally-disturbed 
galaxies uses the mean value of the unmasked areas of this image 
together with a noise correction frame. 




Fig. 5. — Sample galaxy in Abell 795 automatically discarded due 
to extended light profile of bright nearby galaxy. Galaxies whose 
tidal parameters decrease by more than 5% when the mask sizes 
for bright sources are increased by 10% are discarded as this is an 
indication that the high tidal parameter is likely due to extended 
light from neighboring bright sources rather than from tidal fea- 
tures. Left panel: Galaxy image shown with 15" scale bar, center: 
residual image, right: residual with larger masking. Percent change 
in tidal parameter with larger masking: —16%. 

where l Xjy is the pixel value at x,y of the object frame 
and M x ^ y is the pixel value at x,y of the model frame. A 
background-noise correction frame is created by adding 
the background sky image created earlier to the model 
image and then dividing by the model. The resulting 
frame is then median filtered and a tidal parameter is re- 
calculated for the model plus noise frame. The final tidal 
parameter, corrected for noise, is then found as follows: 

T - FF 2 — T 2 (2) 

c y galaxy model V / 

This tidal parameter is particularly sensitive to devi- 
ations when the model values are small. Because our 
models often go to zero toward the image edges, we find 
that we increase the relative stability of our tidal param- 
eters by first dividing the galaxy images by their mean 
background counts so that all images have a pedestal of 
1 count. 

The calculation of a tidal parameter does not by it- 
self identify a set of tidally disturbed galaxies because 
there is no absolute reference. We define a threshold T c 
parameter that is best matched to what we visually iden- 
tify as tidally disturbed galaxies. The correspondence is 
not perfect, as we describe further below, but we settle 
on defining those galaxies with T c > 9 x 10 -4 as tidally- 
disturbed. We find that many galaxies have higher val- 
ues for the tidal parameter because light from the ex- 
tended haloes of nearby bright galaxies extends beyond 
the mask. In an effort to eliminate these false-positives 
we reject from our analysis candidate tidally-disturbed 
galaxies whose tidal parameters decrease by more than 
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Fig. 6. — Sample galaxies flagged as disturbed in the auto-sample, 
but excluded from the cleaned-auto sample in Abell 553 (top) and 
Abell 401 (bottom). Left column: galaxy images displayed with 
15" scale bars, right: region used for the residual calculation that 
has the ellipticity of the primary galaxy and has regions masked. 
Top panel: it is not clear whether these two galaxies are inter- 
acting or are a close superposition. The high mean residual value 
appears to be primarily caused by light from the adjacent galaxy 
spilling beyond its mask rather than any visible tidal structure 
(T = 1.4 x 10 -3 ). Bottom panel: the high residual value appears 
to be primarily due to the extended light profile of a nearby giant 
elliptical galaxy (T = 1.8 x 1CT 3 ). 

5% when the mask sizes for bright sources are increased 
by 10% (e.g. Fig. Ej). While it is true that this filter 
is biased against galaxies in dense regions, it would only 
bias our results if real tidal features are being masked by 
the increased mask sizes. We verified that the fraction of 
discarded galaxies that were visually-disturbed was con- 
sistent (within statistical errors) with the fraction of all 
galaxies that are visually-disturbed within a given radial 
bin, indicating that these filters do not introduce a ra- 
dial trend. The remaining galaxies with T c > 9 x 10 -4 
make up the auto-detected sample of tidally-disturbed 
galaxies. Lastly, because some false positives arise due to 
image quality issues, unmasked sources, and diffuse light 
gradients from nearby bright galaxies, we visually inspect 
this sample to construct our cleaned-auto-detected sam- 
ple of tidally-disturbed galaxies. The procedure we have 
just described consists of many steps, all of which affect 
the number of systems that will ultimately be classified 
as mergers, and thereby complicate comparison to mod- 
els on an absolute scale. However, these filters should be 
mostly insensitive to cluster centric radius, particularly 
outside the crowded central part of a cluster right near 
the brightest cluster galaxy, and as such are not expected 
to affect our results. 

3710 galaxies satisfy all of the automated selection cri- 
teria (cluster color-magnitude band, M r > < —20, r _1//4 
radial brightness profile, successful IRAF bmodel gen- 
eration, and less than 50% of the area surrounding the 
galaxy masked). Of these, we discard 79 that are prob- 
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Fig. 7. — Examples of galaxies marked as disturbed in both the 
visual and cleaned-auto samples in Abell 21 (top) and Abell 1068 
(bottom) . Left column: galaxy images displayed with 15" scale 
bars, right: region used for the residual calculation that has the 
ellipticity of the primary galaxy and has regions masked. Top 
panel: the galaxy is clearly undergoing a merger; note the bright 
shell below the galaxy and the diffuse tidal stream above the galaxy 
(T = 1.2 x 10 -3 ). The mean surface brightness within the red box 
"A" is 25.4 and within box "B" is 25.1. Bottom panel: also a 
clear case of a recent or ongoing merger with a distinct tidal tail 
extending to the lower left of the galaxy and tidal debris to the 
left and above the galaxy (T = 1.1 x 10 -3 ). Galaxies that are 
identified by both the visual and cleaned-auto samples tend to 
have the strongest and most distinctive tidal features. 

lematic due to masking issues and background light gra- 
dients (e.g. Fig. 3]). Of the remaining 3631 galaxies, 202 
were flagged as potentially-disturbed by our tidal param- 
eter threshold, while 109 galaxies were visually-identified 
as tidally-disturbed (see Figures [3 and [8] for examples). 
While there are a large number of disturbed galaxies that 
the automated sample identified that were missed by the 
visual inspection (see Fig. [8] for an example), the sample 
also includes many false identifications (usually due to 
image quality problems or unmasked sources). By eye 
we discard obvious false detections from the automated 
sample to make a cleaned-auto sample of 122 galaxies 
(see Fig. [6j) out of 3551 total galaxies (hereafter referred 
to as "all"). The samples are compared below. 

3. RESULTS 

3.1. Rates of Tidal Features 

The intersection of the visual and cleaned-auto samples 
includes 42 galaxies, while the union of the two samples 
includes 189 galaxies. There is no single characteristic 
that uniquely determines whether a galaxy is identified 
or missed by the visual or cleaned-auto samples. The 
visual sample includes galaxies that have tidal features 
that are close to ellipsoidal, of smaller spatial extent, 
or of such low surface brightness that they are missing 
from the cleaned-auto sample. The cleaned-auto sample 
includes galaxies whose residual images show strong fea- 
tures that were not included in the visual sample because 
it was not evident whether these are actual tidal features, 
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Fig. 8. — Examples of galaxies in either the visual or cleaned-auto 
samples (but not both) in Abell 399 (top and middle panels) and 
Abell 553 (bottom). Left column: galaxy images displayed with 
15" scale bars, right: region used for the residual calculation that 
has the ellipticity of the primary galaxy and has regions masked. 
Top panel: Sample galaxy flagged as disturbed in the visual sample, 
but not in the cleaned-auto sample (T = 8.6 x 10 -4 ). Note the ring 
of tidal debris extending from the top region of the galaxy. The 
mean surface brightness within the red box "A" is 25.7 and within 
box "B" is 26.3. Middle panel: Another sample galaxy flagged as 
disturbed in the visual sample, but not in the cleaned-auto sample 
(T = 7.1 x 1CT 4 ). Note the series of tidal shells. Bottom panel: 
Sample galaxy flagged as disturbed in the cleaned-auto sample, but 
not in the visual sample (T = 9.1 x 10 -4 ). Note the diffuse shell 
extending from the lower left portion of the galaxy. 

spiral arms, or nearby satellite galaxies. There are also a 
few instances of the cleaned-auto sample missing galax- 
ies whose features are very prominent, because these have 
been masked as separate sources by SExtractor. The in- 
tersection sample tends to include the strongest, most 
spectacular features, while the union is a more complete 
sample, but suffers from a larger fraction of questionable 
detections. 

The tidal feature rate among early-type galaxies in 
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these environments, using the cleaned- auto and visual 
samples, is 3.4 ±0.3% and 3.0 ±0.3%, respectively, where 
the uncertainties are purely statistical. The union of the 
two samples results in a rate of 5.3 ± 0.4%. As we men- 
tion in §2.2 a tight color selection could result in a bias if 
interacting galaxies are bluer than the mean. The mean 
color offset from the red sequence is 0.00 ± 0.07 for the 
visually-selected sample and —0.01 ±0.06 for the cleaned- 
auto sample, while the mean color offset for all galaxies 
analyzed is 0.00 ±0.06, demonstrating that our interact- 
ing samples do not suffer from a color bias. The rate of 
tidally disturbed cluster galaxies we find is significantly 
lower than the rate observed by iTal et all (|2009l ). but 
the difference is likely because our sample is ~ an or- 
der of magnitude more distant, resulting in the smaller 
apparent size of galaxies and their associa ted tidal fea- 
tures . Our ra tes match more closel y those of lBridge et al.l 
(|2010T ) and iMiskolczi et al.l (|2011l ). although again it is 
complicated to compare results among studies. 

Although the measured rate of tidal features appears 
small at first, some minor adjustments can lead to an 
inference that the majority of the initial cluster elliptical 
population must have experienced a major merger. First, 
the ~ 3% rate, implies that at least 6% of the initial pop- 
ulation experienced a merger because a merger involves 
two galaxies. If one then considers that the tidal sig- 
natures will not last indefinitely the numb ers rise again . 
If, for example, tidal features last 2 Gyr (jQuinn I fl984) 
then the rate of mergers increases by a factor of about 
5 to ~ 30%. If one further considers that mergers were 
probably more likely in the past, then the number rises 
even more. None of this is empirically well determined, 
but our fundamental conclusion is that our low rate mea- 
surement does not necessarily imply a low incidence of 
mergers over the lifetime of these galaxies. 

3.2. Interlopers and Projection Effects 

Before proceeding to examine the spatial distribution 
of galaxies with tidal features, we stop to consider the 
effect of interloping galaxies on our radial trends. With- 
out spectroscopic redshifts, we do not know with cer- 
tainty whether or not a particular galaxy is in the clus- 
ter. This acts to muddle any real correlations with ra- 
dius. To estimate the effect of interlopers within our 
color-selected samples, we count the number of galaxies 
selected using color-magnitude bands of the same size 
displaced redward by 0.26 mag (the thickness of our se- 
lection area). These displaced color- magnitude bands 
contain 2834 galaxies, while the bands centered on the 
red sequences contain 11,904 galaxies, suggesting that 
approximately 80% of our selected galaxies are cluster 
members. We also count the n umber of gal axies in the 
four deep fields of the CFHTLS (jGwvn 112008ft that would 
satisfy our galaxy selection criteria. The count from the 
CFHTLS-Deep fields yielded a much lower estimation of 
the number of interloping galaxies, but had large field- 
to-field variation. Therefore, we conservatively use the 
interloper estimate from the displaced color-magnitude 
bands and repeat our tidal parameter analysis for these 
redder galaxies to determine what effect they may have 
on the results. A higher fraction of the background galax- 
ies are ellipticals, because we are dealing with brighter 
and redder galaxies, resulting in a slightly higher frac- 
tion of these (34% versus 30% of the cluster sample) 



meeting all selection criteria and contributing to the fi- 
nal results. As expected for a background population, 
the fraction of interlopers that are tidally disturbed does 
not correlate with clustercentric radius — demonstrat- 
ing that crowding toward the center of the cluster does 
not artificially create a radial effect. However, interlop- 
ers are a larger fraction of our "cluster" sample at larger 
clustercentric radii (as high as 45% beyond R200 vs less 
than 10% within R200) due to the steep decrease in the 
number density of actual cluster galaxies at larger clus- 
tercentric radii. The background galaxies have slightly 
higher rates of tidal features (4. 3^0*. 7% visually-detected 

and 4.5^0.7% auto-detected, with the statistical uncer- 
taint ies calcu l ated using the Poisson single-sided limits 
from Gehrels (1986)) than the cluster sample (3.1 ±0.3% 
visually-detected and 3.3 ±0.3% auto- detected) suggest- 
ing that they might slightly influence any radial behavior. 

3.3. Clustercentric Radius 

To examine the rates of tidal features as a function of 
projected clustercentric radius, we rescale th e projected 
radii u sing the R200 value of each cluster from lSand et al.l 
(2012). These radii were estimate d using the conversion 
betwe en Lx and M200 found by iReiprich fc Bohringerl 
(2002) and the correspondence between M200 and i?200- 
Clustercentric radii are measured from the brightest clus- 
ter galaxy (though we verified that the results do not 
change significantly when the radii are measured from 
the X-ray center). 

Not all clusters are sampled to the same normalized 
projected clustercentric radius, so we consider different 
subsamples. We create a first subsample of the 1657 
galaxies within R200 in the 50 clusters for which our sur- 
vey is complete out to i?200, a second subsample of the 
1540 galaxies within 2^200 m the 26 clusters for which 
we are complete to 2i?200, and a combined sample of all 
galaxies within R200 in clusters complete to R200 and all 
galaxies within 2i?200 in clusters complete to 2i?200- In 
Figure [9] we plot the relative radial distributions of "all" 
galaxies and our various samples of galaxies with tidal 
signatures. 

The distributions shown in Figure [9] make clear that 
there is a significant deficit of galaxies with tidal signa- 
tures towards the centers of these clusters. The deficit is 
evident within ~ 0.5i?200, which is well outside an area 
about the brightest central galaxy where one might ex- 
pect that confusion could lead to lower detection rates. 
In fact, as we find with a control sample of background 
galaxies that shows no gradient in tidal features, we are 
not affected by confusion in the center of these clusters 
(§3.2). Finally, these findings are independent of whether 
visual or automated techniques are utilized. 

To place these conclusions on a more solid statistical 
footing, we apply K-S tests to evaluate the likelihood 
that the cumulative distributions of galaxies with tidal 
features are drawn from the parent distribution of "all" 
galaxies. For the sample complete within R200 we find 
that the tidal sample was not drawn from the parent 
sample with 98.2% confidence for the visually-selected 
sample, 91.5% confidence for the automated selection, 
94.5% for the intersection of the two samples, and 97.4% 
for the union. For the sample complete to 2^200, the 
confidence limits fall below 2<j significance. 
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Fig. 9. — Histograms of galaxies versus projected clustercentric radius. The left section is the sample of 50 clusters for which the survey 
is complete within -R200, whereas the right section is the sample of 26 clusters complete within 2i?200- Top panel: "all" galaxies, 2nd panel 
down: visually-selected disturbed galaxies, 3rd panel down: cleaned-auto sample of disturbed galaxies, 4th panel down: union of visually 
and cleaned-auto-selected galaxies, bottom panel: intersection of visually and cleaned-auto-selected galaxies. 



The radial variations in the incidence of tidal signa- 
tures suggests that either (or both) the creation or de- 
struction of tidal features depends on clustercentric ra- 
dius. More appropriately, the variation depends on a 
physical characteristic that drives the phenomenon and 
which is itself correlated with clustercentric radius. A 
straightforward hypothesis is that tidal signatures are 
more easily erased in a dense environment due to shorter 
dynamical times. While this is likely to be true near the 
center of the clusters, where the intracluster light is built 
up, by the time one is considering galaxies at R200 and 
beyond, it may seem unlikely that the global potential 
will strip tidal material. 

To tes t this hypothe s is, we rely on numerical simula- 
tions by Ru dick et al.l (|2QQ9[ ) that find that the decay 
time of tidal streams is approximately 1.5 times the dy- 
namical time at their location in the cluster. The dy- 
namical time at clustercentric radius r is given by, 



tdyn 



GM 



(3) 



where M is the mass enclosed within r. We therefore set 
the lifetime of tidal features, /(r), to be l.btd yn - We de- 
termine the enclosed mass usin g the NFW density profile 



t ermine tn e enclosed mass using t 
(N avarro. Frenk. fc Wh ite 1997): 



p(r) = 
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(r/a)(l + r/a) 2 



(4) 



where a = R200/C and we take c = 2.6 to be an appropri- 
ate concentrati on parameter for the di stribution of galax- 
ies in clusters (jBudzynski et al.ll2Q12f ). If the fraction of 
galaxies to form tidal features, //, is uniform with clus- 
tercentric radius and time, then the expected fraction of 
galaxies observed to have such features would be: 
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f f !$ for l(r) < t 
f f for l(r) > t , 



(5) 



where to is the age of the Universe. The observed fraction 
as a function of projected clustercentric radius is then 
given by calculating the average fraction along the line 
of sight: 



/06s {^proj ) 



J-oo fobs(x)n(x)dx 
JZo n{x)dx 



(6) 



where x = 



— Tp ro j and n is the number density of 

our "all" sample. Because the NFW profile is not a good 
description of the mass or galaxy number density to ar- 
bitrarily large radius, we replace the integration limit 
with ±10i?200- Using this prescription and assuming 
that galaxies trace mass, we find that the distribution 
of galaxies predicted as a function of projected radius by 
the NFW profile roughly reproduces our observed distri- 
bution (see Fig. [TO]) . 

We now compare our results to the expected distri- 
bution of observed tidal features if the lifetime of the 
features is proportional to the local dynamical time and 
the generation of tidal features is independent of time 
or clustercentric radius (see Fig. [TT]) . While the cumu- 
lative distribution of tidally-disturbed galaxies still lies 
below the model prediction within i?200, the deviation is 
still only statistically significant for the visually-selected 
sample of disturbed galaxies. For the sample complete 
within i?200> the K-S test shows that the visually-selected 
sample is drawn from a different distribution than the 
"all" sample with 95.8% confidence, the automated se- 
lection differs with 81.3% confidence, the intersection of 
the two samples is discrepant at 91.0% confidence, and 
the union is discrepant at 91.7%. For the sample com- 
plete to 2i?200 none of the samples of disturbed galax- 
ies are discrepant from the "all" sample at the 2a level. 
These results are robust with reasonable changes in the 
concentration parameter used. If we instead use c=6, 
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Projected Clustercentric Radius (R/R 200 ) 

Fig. 10. — Cumulative fraction of observed galaxies (within clus- 
ters with FOV complete to 2i?20o) as a function of projected clus- 
tercentric radius (solid blue line) compared to cumulative distri- 
bution of galaxies modeled by an NFW profile with c=2.6 and in- 
tegrated out to a maximum radius of 10i?200 (dashed green line). 



which is more accura te for the dark matter distribution 
for ha los of this mass (Ma ccio. Dutton. fc van den Boschl 
2008), the K-S statistics are largely unchanged. Addi- 
tionally, we tested a range of coefficients (0.5 — 2Md yn ) 
for the decay time behavior of tidal streams and, based on 
the one to two percent changes on the K-S probabilities, 
affirm that our conclusions are independent of modest 
changes to the choice of the coefficient. We find that this 
toy model, in which the lifetime of the features is pro- 
portional to the local dynamical time and the generation 
of tidal features is independent of time or clustercentric 
radius, qualitatively agrees with the observed deficit of 
tidal features at small clustercentric radii, but possibly 
would not solely account for the magnitude of the deficit 
within ~ 0.5i?200- 

3.4. Local Density 

When considering the effect of environment on clus- 
ter galaxies there is always a tension between the role 
of global vs. local environment. Part of the difficulty in 
addressing such issues is that the two measures of envi- 
ronment are correlated, and therefore large samples are 
needed to even begin to attempt to disentangle the two. 
A second difficulty is that local environment is inher- 
ently difficult to measure in clusters due to the compli- 
cated line-of-sight structure and relatively poor distance 
estimates. Here we make an attempt to study the de- 
pendence of tidal features on local density, but we face 
these same difficulties. 

We quantify the local density using the standard mea- 
sure of the projected distance to the 10th nearest galaxy 
with M r < — 18 that is also within the cluster's red se- 
quence in an attempt to mitigate contamination. We find 
no evidence that the rate of tidal features is dependent 
on local density as both the visual and the cleaned-auto 
samples differ from the full cluster sample at less than 2a 
significance using a K-S test. We next attempt to remove 
the correlation between local density and clustercentric 
radius by dividing the galaxies into clustercentric annuli 
and within each annulus ranking galaxies uniformly by 
local density to range from 1 to 10 (10 being most dense). 
We then sum the density ranks of disturbed galaxies for 



all distance bins and compare that distribution to a ran- 
dom distribution (see Fig. [T2]) . By doing this, we are 
asking whether at each radius the galaxies with tidal fea- 
tures are preferentially in environments of either low or 
high local density. We find that disturbed galaxies have 
no statistically significant deviations in ranks from a ran- 
dom distribution. 

Although we find no evidence for local density having 
an effect on the incidence of tidal features, we caution 
that this is simply an absence of evidence rather than 
evidence of an absence, particularly because our mea- 
surement of local density is so crude. 

4. CONCLUSIONS 

In a sample of 54 galaxy clusters (0.04 < z < 0.15) 
containing 3551 early- type galaxies suitable for study, 
we identify those with tidal features both interactively 
and automatically. This constitutes the largest sample 
studied to date for signs of environmental dependence in 
the incidence of tidal features. We find tidal features in 
~ 3% of galaxies in this sample, with data of this par- 
ticular depth. Regardless of the method used to classify 
tidal features, or the fidelity imposed on such classifica- 
tions, we find a deficit of tidally disturbed galaxies with 
decreasing clustercentric radius that is most pronounced 
inside of ~ 0.5i?200- 

Although this trend could be attributed to a rise in 
galaxy-galaxy interactions at large clustercentric dis- 
tances, where the galaxy pair velocities might be better 
tuned to mergers, an alternative interpretation is that 
tidal features are preferentially erased at small cluster- 
centric radii. We examine this hypothesis with a simple 
toy model that links tidal feature survival time to the lo- 
cal dynamical time and find that although qualitatively 
the effect is in the correct sense to explain the data, the 
model falls somewhat short quantitatively. Nevertheless, 
given the relative success and the limited statistics, we 
cannot exclude this model as a possible explanation for 
the radial trend in the incidence of tidal tails. 

We also search for a dependence of the incidence of 
tidal features on local density and find no statistically 
significant evidence for such, even after accounting for 
the correlation between clustercentric radius and local 
density. Unfortunately, the large uncertainties in the 
measure of local environment limit the interpretation of 
a null result. 

We demonstrate that interesting behavior exists in the 
rate of tidal features among cluster early-types as a func- 
tion of clustercentric radius. This measure should pro- 
vide some guidance for models of merging among early- 
types, which is conjecture d to be sign i ficant in the build- 
up o f the red sequence (Be ll et all 120041 : iFaber et~aT1 
2007). Models of merging in and around clusters should 
aim to provide merger estimates based on observable fea- 
tures rather than theoretical constructs. We expect that 
such work will provide strong constraints on the effect of 
the cluster environment on the structure of galaxy ha- 
los, the build-up of the red sequence of galaxies, and the 
origin of the intracluster stellar population. 
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Fig. 11. — Cumulative fraction of galaxies versus projected clustercentric radius. The left panel is the sample complete within R200, 
the center panel is the sample complete within 2i?200, and the right panel is the combination of the two samples (all galaxies within R200 
in clusters complete to R200 and all galaxies within 2i?200 in clusters complete to 2i?20o)- Solid black line: "all" galaxies (dominated by 
undisturbed galaxies), dotted purple line: cleaned-auto sample, short-dashed red line: visually selected galaxies, long-dashed cyan line: 
union sample, dot-dashed green line: intersection sample, solid orange line: expected distribution for model (Eqn. [6} in which lifetime 
of the features is proportional to their dynamical time with a uniform incidence of features at all radii and an overall fraction of tidally 
disturbed galaxies set to match the observed value. All three panels show a deficit of tidally disturbed galaxies within ^200 • 




2 4 6 8 10 

Ranked Density Bin 



Fig. 12. — Histograms of galaxies versus ranked local density in 
radius annulus (a random distribution would be flat across all den- 
sity bins). Top panel: visual sample of disturbed galaxies, 2nd 
panel down: cleaned-auto sample, 3rd panel down: union of visu- 
ally and cleaned-auto-selected galaxies, bottom panel: intersection 
of visually and cleaned-auto-selected galaxies. 



on observations obtained with MegaPrime/MegaCam, 
a joint project of CFHT and CEA/DAPNIA, at the 
Canada- France-Hawaii Telescope (CFHT) which is oper- 
ated by the National Research Council (NRC) of Canada, 
the Institute National des Sciences de l'Univers of the 
Centre National de la Recherche Scientifique of France, 
and the University of Hawaii. This work is based in 
part on data products produced at TERAPIX and the 
Canadian Astronomy Data Centre as part of the Canada- 
France-Hawaii Telescope Legacy Survey, a collaborative 
project of NRC and CNRS. 
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